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ABSTRACT 

Quasi-periodic signals have yielded important constraints on the masses of 
black holes in galactic X-ray binaries, and here we extend this to active galactic 
nuclei (AGN). We employ a wavelet technique to analyze 19 observations of 
10 AGN obtained with the XMM-Newton EPIC-PN camera. We report the 
detection of a candidate 3.3 kilosecond quasi-period in 3C 273. If this period 
represents an orbital timescale originating near a last stable orbit of 3 Rs, it 
implies a central black hole mass of 7.3 x 10^ Mq. For a maximally rotating black 
hole with a last stable orbit of 0.6 Rs, a central black hole mass of 8.1 x 10^ 
Mq is implied. Both of these estimates are substantially lower than previous 
reverberation mapping results which place the central black hole mass of 3C 273 
at about 2.35 x 10^ Mq. Assuming that this reverberation mass is correct, the 
X-ray quasi-period would be caused by a higher order oscillatory mode of the 
accretion disk. 



Subject headings: galaxies: individual(3C 273) 
galaxies 



galaxies: active — X-rays: 



Introduction 



Quasi-periodic oscillations (QPOs) are thought to originate i n the inner accreti on disk 
of a black hole or neutron star in an X-ray binary (XRB) system ( Ivan der Klisll2000l ). Con- 
sequently, QPOs have been used in galactic XRBs to introduce important constraints on the 
masses of the central black holes of these systems. 

Previous work has revealed that AGN and XRBs are alike: noise power spectra have 
shown that similar phys i cal processes may be underlying the X-r a y variability in both 
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McHardy et al.l |2004| . |2005| ) . Taking this resemblance into account and assuming that ac- 
cretion onto a stellar-mass black hole is comparable to accretion onto a supermassive black 
hole, one would expect some AGN to exhibit QPOs similar to those observed in XRBs. In 
supermassive black holes (10^-10^ Mq), these QPOs would be at much lower frequencies 
than those we find in stellar-mass black holes (~ 10 Mq). Low frequency quasi-periods (LF 
QPOs) in XRBs range from 50 mHz to 30 Hz; scaling from a ~ 1 H z QPO in a 10 M^r^ XRB, a 
LF QPO in an AGN would occur on timescales of days to months (IVaughan fc Uttleyll2005l ). 
too long to be detectable for the AGN in our sample. On the other hand, high frequency 
QPOs (HF QPOs) in XRBs have values of > 100 Hz and assuming a 1 / Mb h scaling of fre- 
quencies, fuFQPO ~ 3 X 1Q^^{Mbh I^^^Mq)"^ Hz (lAbramowicz et al.ll2004l ). corresponding 
to timescales greater than 400s for AGN. 

While this parallel between AGN and XRBs seems promi sing, no claim of a n X-ray 
quasi-period in an AGN has been found to be statistically robust. IVaughan fc Uttleyl (l2005bl ) 
remark that a major source of false detections arise from assuming an inappropriate back- 
ground noise power spectrum. X-ray variations of AGN have intrinsically red noise power 
spectra (i.e. the power spectra have a continuum resembling a power law with a steep slope; 
Press 1978), however many purported QPOs in AGN are compared against an assumed 
background of white noise (i.e. Poisson photon noise or a fiat spectrum). For example, 
in a ~5 day ASCA observation of IRAS 18325-5926 the sign ificance of the candidate peri- 
odicity was estimated w ith white noise (jlwasawa et al.lll998l ). After including red noise in 
the periodogram fitting. I Vaugharu (|2005[ ) found that the candidate periodicity was no longer 



significant at the 95% level. iFiore et al.l (119891 ) also claimed high (> 99%) significance peaks 
in NGC 4 151, however, after fitting red noise and Poisson photon noise components of the 
spectrum IVaughan fc Uttleyl (l2005bl ) showed that the significances of the QPOs fell below 
the 95% confidence level. It is also difficult to constra in the significance of possible QPOs 
due to power spectra effects JVaughan fc Uttlevll2005bh. EXOS AT data of NGC 5548 were 



report ed to have a significant period (jPapadakis fc Lawrencd Il993l ). but iTagliaferri et al. 



(119961 ) later showed that the significance of the candidate QPO was lower than previously 
reported once the uncertainties in modeling the spectrum were taken into consideration. 

This lack of statistically significant evidence for QPOs in AGN has led to questions of 
whether existing X-ray observations of A GN are sensitive enough to detect QPOs even if 
they are present (IVaughan fc Uttleyl l2005l ). 



Here we use a different technique to search for significant periodic structures in the time 
variability data that have been collected for AGNs with XMM-Newton. We use a wavelet 
transform technique, which can have certain advantages relative to periodograms and Fourier 
power spectra, the methods that have previously dominated the literature. The wavelet 
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technique, which has become widely used in other branches of science, is particularly useful 
in identifying signals where the period or its amplitude changes with time. This technique is 
applied to the XMM-Newton data from bright 10 AGNs, with special care taken to properly 
treat the noise characteristics and error analysis, and we find a candidate 3.3 ks quasi-period 
in 3C 273. 

In Section 2, we present our observations and data reduction steps. In Section 3, we 
provide an overview of the two wavelet techniques used in our analysis: the continuous 
wavelet transform and the cross-wavelet transform. The results of these two techniques as 
well as significance tests are presented. We also discuss structure function analysis for the 
AGN in our sample. In Section 4 we argue that this 3.3 ks quasi-period in 3C 273 is consistent 
with what we would expect from oscillations in the accretion disk around the supermassive 
black hole based on current black hole mass estimates. 



2. Observations and Data Reduction 

The 10 AGN in our sample were selected because they are bright and have XMM- 
Newton EPIC-PN camera observations which exceed 30 kiloseconds (ks). In total, we have 
19 observations and each observation's ID, date, length, and average counts are listed in 
Tabled! All observations are in the energy range 0.75 to 10 keV and most were taken in 
small window mode, which has a readout time of 6 milliseconds (ms). The only exception is 
NGC 4151 Observation ID (Obs. ID): 0112830201, which was taken in full frame mode with 
a readout time of 73.4 ms. 

Observation Data Files (ODFs) were obtained from the on-line XMM-Newton Science 
Archive and later reduced with the XMM-Newton Science Analysis Software (SAS, v. 7.0.0, 
6.1.0, 5.4.1). Source light curves, with 5 s bins, were extracted for a circular region centered 
on the source (~20"). Background light curves were obtained from a nearby rectangular 
source-free region and subtracted from the source light curves. These rectangular background 
regions were larger than the source regions and were accordingly scaled down. Due to strong 
flaring, the last few kiloseconds of data are excluded from most observations. The count 
rates for the target sources are orders of magnitude greater than the background count rate 
in the detection cell, so a rise in the background is unimportant. We removed these last 
few kiloseconds of data from the data stream just to be very cautious. We note that in the 
observation of 3C 273 with the claimed detection, including the periods with flaring does 
not change our results. 

Some of the observations in our sample are affected by pile-up. Pile-up occurs when more 
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than one X-ray photon arrives in a pixel before the pixel is read out by the CCD, making 
it difficult to distinguish one high energy photon from two lower energy photons. Pile-up 
can also occur when photons striking adjacent pixels are confused with a single photon that 
deposits charge in more than one pixel. Depending on how many pixels are involved, this 
is called a single-, double-, triple-, or quadruple- pixel event. The SAS task EPATPLOT 
measures the pile-up in an observation and the results for our target with the highest count 
rate, MKN 421 (Table 1), are shown in Figure [H When we compare the expected fractions 
of pixel events (solid lines) with those actually measured in the data (histograms) for the 
range 0.75 to 10 keV we see that a larger than expected fraction of double events (third 
histogram from top, dark blue in electronic edition) is measured as well as a larger fraction 
of triple and quadruple events (bottom two histograms), although to a lesser degree, while 
single events (second histogram from top) are lower than expected, indicating the presence 
of pile- up. Pile-up leads to a general reduction in the mean count rate as well as a reduction 
in the magnitude of variations. We will explore the influence of pile-up on our data in more 
detail when we discuss structure functions in Section 13.2. 1[ 



3. Data Analysis and Results 
3.1. Wavelet Analysis 



3.1.1. The Continuous Wavelet Transform 



The continuous wavelet transform (CWT) is the inner product of a dilated and trans- 
lated mother wavelet and a time series f{t), the idea being that the wavelet is applied as a 
band-pass filter to the time-series. The continuous wavelet transform maps the power of a 
particular frequency (i.e. dilation) at different times in translation-dilation space, giving an 
expansion of the signal in both time and frequency. Hence, the continuous wavelet transform 
not only tells us which frequencies exist in the signal, but also when they exist, allowing us 
to see whether a timescale varies in time. This is the wavelet technique's advantage over 
Fourier transforms in detecting quasi-periods. In addition, the Fourier transform is not 
suited for detecting quasi-periods since non-periodic outbursts will spread power across the 
spectrum and windowing will cause power to appear at low frequencies, potentially obscuring 
quasi-periodic signals. 

Throughout this paper, we follo w iHughes et al. Jl998h and iKellv et al.l J2003h and ref- 



erences within. In previous studies ( Hughes et al. 19981 : Kelly et al. 



2003 



Liu et al. 2005 
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Kadler et al.ll2006l ) we have found the Morlet wavelet 

'ipMorlet = 



l/4gifc^tg-|i2|/2 



with = 6 to be an excellent choice. The value of is a satisfactory compromise between 
a value small-enough that we have good resolution of temporal structures, and large-e nough 
that the admissibility condition is satisfied, at least to machine accuracy (lFargdll992l ). The 
wavelet, being continuous and complex, permits a rendering in transform space that high- 
lights temporally localized, periodic activity - oscillatory behavior in the real part and a 
smooth distribution of power in the modulus - and being progressive (zero power at nega- 
tive frequency), is optimal for the study of ca usal signals. W e have deliberately avoided any 
form of weighti ng, such as that introduced by iFosten (119961 ) to allow for uneven sampling, or 
Johnson! (120061 ) to rescale within the cone of influence, in order to facilitate our interpretation 
of the cross wavelet, and to allow the use of existing methods of significance analysis. 

From this mother wavelet, we generate a set of translated {t') and dilated (/) wavelets 



V',,,(t) = i=^(^),ZG3ft+,te3fJ 



(2) 



and we then take the inner product with the signal F[t) to obtain the wavelet coefficients 



fm!At)dt. 



(3) 



The wavelet coefficients are later mapped in wavelet space which has as coordinates trans- 
lation and dilation, and so periodic behavior shows up as a pattern over all translations at 
a specific dilation. 

By way of example. Figure [2] shows the real part and the power of the continuous wavelet 
transform (second and bottom panel, respectively) for a sinusoidal signal of varying frequency 
(top panel). Here, the real part of the transform shows oscillatory behavior corresponding 
to the two periodicities of the sinusoidal signal at dilations of 3s and 6s with a break in 
translation at 50s corresponding to the time where the change in frequency occurs. The 
bottom panel in Figure [2] shows that the power of the continuous wavelet transform is 
concentrated at these two frequencies as well. 

The hatched area in both panels of Figure [2] represents the cone of influence: the region 
where edge effects become important. It arises because discontinuities at the beginning and 
end of a flnite time series result in a decrease in the wavelet coefficient power. Also shown 
in the header of Figure [2] are the number of dilations used (A^^^) and the ranges of dilations 
explored. We discuss a and the normalization of Figure [2] in the following section. 
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3.1.2. Significance Tests 
gnificance tests can be cr eated for the continuous wavelet transform and here we 



Torrence fc Compd (119981 ) . First, one compares the wavelet power with that of an 



appropriate background spectrum. 

We use the univariate lag-1 auto-regressive [AR(1)] process given by 

where a is the assumed lag-1 autocorrelation and Zn is a random deviate taken from white 
noise. Note that a = gives a white noise process. Throughout this paper, we will use 
"white noise" to refer to an AR(1) process with a = 0. Red noise is sometimes used to refer 
to noise with a = 1, however, throughout this paper we apply the term to any non-zero a. 

The normalized discrete Fourier power spectrum of this process is 

p. = f5) 

^ 1 + -2acos{2'KSt/Tj) ^ ' 

where Tj is the associated Fourier period for a scale Ij. We use the above two equations to 
model a white noise or red noise spectrum. 

The global wavelet power spectrum (GWPS) is obtained by averaging in time 

flil^ = ^>i=^^fil^.9?■ (6) 

Here, ij and i'^ are the indices of the initial and final translations t[ outside of the cone 
of influence at a given scale Ij. Nj is the number of translations t'^ outside the cone of 
influence at that scale. Assuming a background spectrum given by Eqn. [5] we estimate the 
autocorrelation coefficient (a) by calculating the lag-1 and lag-2 autocorrelations, ai and 
a2. The autocorrelation coefficient is then estimated as a = (ai + ^foq)/'!. The background 
spectrum Pj then allows us to compute the confidence levels. It is assumed that the time 
series has a mean power spectrum given by Eqn. [5] and so if a peak in the wavelet power 
spectrum is significantly above this background spectrum, then the peak can be assumed to 
be a true feature. If the values in the time series fit) are normally distributed, we expect 
the wavelet power |/p to be distributed with two degrees of freedom (xD- The square 
of a normally distributed variable is distributed with one degree of freedom and the 
second degree of freedom comes from the fact that both the real and imaginary parts of the 
complex / are normally distributed. For example, to determine the 95% confidence level, 
one multiplies the background spectrum (Eqn. [S]) by the 95th percentile value for x^. In 
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Figure [3] we show the GWPS of a time series along with 99% and 95% confidence level for 
a red noise process and 99% confidence level for a white noise process. 

The distribution for the local wavelet power spectrum is 

where the arrow means "distributed as," is the variance, and v is the number of degrees of 
freedom, which is two here. The indices on the scale / are j=l,2,...,J where J is the number 
of scales, and the indices on the translation t' are i=l,2,...,Niiata- We evaluate this equation 
at each scale to get 95% confidence contour lines and in this paper our continuous transforms 
are normalized to the 95% confidence level for the corresponding red noise process. Doing 
this allows one to see the strength of the wavelet coefficients relative to the 95% confidence 
level of a red noise process. 



3.1.3. The Cross-Wavelet Transform 

Although the continuous wavelet transform is useful in examining how a time series 
varies in time and scale, it does not tell us how the time series varies in dilation over a 
range of scales when assigning a characteristic timescale. Since a quasi-periodic signal has 
no unique dilation we use the cross-wavelet transform (XWT) whi ch filters ou t noise and 



reveals the QPO more clearly. Here we use the XWT introduced by lKelly et al.l (120031 ). 



After the continuous transform identifies that a periodic pattern exists in the data, the 
dilation that characterizes this period is obtained from the global wavelet power spectrum 
and is used to create a sinusoidal mock signal. The continuous wavelet transform coefficients 
of the data signal fa{t) are then multiplied by the complex conjugate of the continuous 
transform coefficients of a mock signal fmif)- The results are mapped out in wavelet space 
and analyzed for a correlation. 

The cross-wavelet transform takes the form 

Ul,t') = JaMr^M (8) 

where the continuous wavelet coefficients fa and fm are given by Equation [31 

Figure H] shows the cross- wavelet for the same sinusoidal signal of varying frequency used 
in Figure [21 The mock signal was calculated using the 6s period found in the wavelet power 
spectrum (see Fig. [3]) and as the concentrations in the real and power panels of Figure [H 
show, the cross-wavelet finds that this 6s period exists in the first half of the time series, 
illustrating the cross- wavelet's ability to highlight a QPO. 



The reader may refer to 
used here. 



Kelly et al, 



(I2OO3I ) for a full review of the cross-wavelet technique 



3.2. Structure Function Analysis 

Since the global wavelet power spectrum compares the observed signal to expected levels 
of red noise and white noise, we created structure functions (SFs) for each of our observations 
to see which noise process dominates the signal at different times. 

A structure function calculates the mean deviation of data points, providing an alternate 
method of quantifying time variations. Here we use a structure function of the first-order 



( ISimonetti et al.lll985l ): 

SF{5t) =< [F{t) - F{t + 5t)f > (9) 

where F{t) is the flux at time t and 6t is a time lag. The slope a of the SF curve in 
log{SF) — log{5t) space depends on the noise processes underlying the signal, giving us an 
indication of the nature of the process of variation. If a = 1 red noise dominates, and for 
flatter slopes of a = Poisson photon noise is significant. A plateau at short time lag is due 
to measurement noise. The transition from plateau to power-law in the structure function 
curve determines where the dominant underlying noise process changes in the object. The 
point of turnover from power-law to plateau at longer time lags corresponds to a maximum 
characteristic timescale. 



3.2.1. Effects of Pileup 

We measure the presence of pile-up in our observations by using the SAS task EPAT- 
PLOT and find that the majority of our sources show varying degrees of pile-up. For example, 
as previously shown in Section [2], MKN 421 Obs. ID: 0099280101 has a modest amount of 
pile-up (see Figured]). In the structure function of this observation (left panel. Figure [5]), 
the flat portion of the structure function curve should have a value of log{SF) = 1, which 
corresponds to the Poisson photon noise inherent in the photon statistics. However, here it 
falls below the Poisson photon noise level. To remove the pile-up we exclude the central core 
of the source in the event file since pile-up is more likely to occur here. For this subtracted 
data, the EPATPLOT output indicates that there is no pile-up and the SF curve is then at 
the expected value for Poisson photon noise (right panel. Figure [5]). Pile-up affects the SF 
because it lowers the overall count rate and thereby Poisson photon noise is underreported. 

We correct for pile-up in the rest of our data by adding a fixed value to log{SF), moving 
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the flat part of the structure function curve up to 1. All of our observations had less than 5% 
pile-up except for PKS 2155-304 Obs. ID 124930301 (6.5%) and both observations of MKN 
421 (~ 10%). Overall, the percentage of pile-up in our sample increases with the number of 
counts except for NGC 4151 Obs. ID 112830201 which is 5% piled-up and has an average of 
only 25 counts. 

3.3. Results 

3.3.1. Wavelet Analysis Results 

Of the observations that we analyzed, only one showed a quasi-period of interest (at 3.3 
ksec), and this occurred in an observation of 3C 273 (ID 126700301). The continuous wavelet 
transform result for this observation is shown in Figure [H] with the quasi-period circled in 
the real and power plots (second and third panel, respectively). One can see that the quasi- 
period appears in the last two-thirds of the observation. In the real plot, the concentrations 
match up with peaks in the light curve, and the power is concentrated at 4.2 x 10^ s. The 
wavelet is sampled with 220 dilations (A*";) ranging between ~ 207.2 s and 2.3 x 10"^ s. We 
note that the data in Figure O are binned from 5 s to 100 s for clarity and that we only 
show the first 56 ks due to background fiaring at the end of the observation. We note that 
including the periods with background fiaring does not change our results. The a found 
from autocorrelation analysis for the unbinned data is 0.14 and this value is used to reach 
the conclusions in this paper. 

The 3.3 ks quasi-period is also evident in the Global Wavelet Power Spectrum (GWPS, 
Figure [7]), which is calculated by summing up the wavelet power spectra at all times. In 
searching for quasi-periodic behavior we excluded time scales above 25% of the time series 
length, where, using spectral methods, too few periods to provide a convincing result would 
be present, and where the cone of infiuence becomes important for the wavelet coefficients. 
On short time scales, experience has shown that sources often exhibit a broad distribution 
of power, with local maxima not well-separated from the mean power level. We selected 
a lower bound for our search, by visual identification of such behavior in the GWPS, in 
conjunction with a concomitant change in behavior of the SF. The solid line in Figure [7] is 
the power spectrum of the signal, which is compared to the power spectrum of white and 
red noise random processes (broken lines). One can see that the 3.3 ks detection exceeds 
the expected levels of white and red noise at the 99% significance level, i.e. the probability 
of the detection is higher than 99% of the noise random processes (the significance of this 
signal is 99.979% relative to red noise with a = 0.14). The origins of the white and red noise 
power spectra were discussed in Section 13.1.21 
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The cross-wavelet analysis for 3C 273 (Fig. [8]) supports the conclusion that a period 
of 3.3 ks is indeed present. Here, the XWT (see Section [3.1.3P compares a mock sinusoidal 
signal with a period of 3282 s with the 3C 273 light curve. The concentration in the cross- 
wavelet transform shows that the 3.3 ks signal is present throughout the observation. As one 
can see, by comparing the crosswavelet signals in juxtaposed bands, the 3.3 ks periodicity 
can be traced over the entire interval. In the CWT (Figure [6]) the 3.3 ks signal is particularly 
strong at late times, and so, due to the limited dynamic range of the rendering, is not evident 
early in the time interval in that figure. 

This periodicity is not detected in the other three observations of this object. In the 
58 ks (Obs. ID 159960101) and 60 ks (Obs. ID 126700801) observations of 3C 273, there 
is a signal at 5000s, but it does not rise above the 99% red noise confidence level (Fig. 
We note that a Fourier analysis of 3C 273 yielded a feature at 3.3 ks, but with a lower 
significance (< 3a) than is found with the wavelet technique. 

We performed Monte Carlo simulations in order to estimate the probability that the 
wavelet technique would claim a spurious detection. As a baseline, we created one thousand 
simulated light curves for Poisson photon noise (Fig. 11) to represent random observational 
errors i.e. photon counting statistics. The simulated light curves were 56 ks long with 5 s 
intervals and we multiplied the mean deviate z„ by 40 to produce an average spread in the 
y-axis of 40 counts to resemble the 3C 273 light curve. Most of the false detections occur 
at timescales less than 2000s which corresponds to 3.6 % the length of the observation and 
supports our earlier point that one can select the lower limit to search for periodicities by 
visual identification of broad distributions of power on short time scales in the GWPS. On 
average, the wavelet technique claims a detection (at or above the significance level reported 
by the wavelet analysis for 3C 273) 0.4% of the time (Fig. [TT]) . 

The Monte Carlo simulations suggest a significantly higher rate of false detections than 
is implied by the statistics based on the GWPS. However, they are consistent with the 
latter estimates within the margin of error, given that only 1000 realizations of a time series 
were generated. Better simulation statistics could be achieved by increasing the number 
of time series realizations by several orders of magnitude, but devoting time and resources 
to this is not warranted. Visual inspection of the simulated light curves reveals that they 
differ qualitatively from the actual time series: a better correspondence can be achieved 
with the addition of randomly distributed Gaussian-profile bursts of fixed, small amplitude. 
Evidently, the process under study is not strictly a stationary, first order one, and the 
formal statistical measures of significance should be regarded as only indicative of the high 
likelihood of a quasi-periodic phenomenon in this source. A more detailed analysis, allowing 
for nonstationary processes, is beyond the scope of this paper. While we have performed 19 
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independent experiments and found only 1 detection we point out that of our 19 data sets 
only 7 have average counts (Table 1) equal to or more than the observation in which we find 
the QPO. One cannot expect to see with equal likelihood, a periodicity of equal strength in 
these weaker AGN. 

We note that independently, the XWT finds evidence for power throughout the obser- 
vation at 3.3 ks (Fig. 8). We measured the 3.3 ks signal strength across the time series from 
the power plot of Figure 8. The power of the 3.3 ks signal is ~ 4000 times stronger than 
shorter and longer dilations, illustrating that the 3.3 ks period is well-constrained. We also 
ran the XWT on this time series with analyzing signals of 2.3 ks and 4.3 ks. The average 
power of these signals is ~ 2 times less than the average power of the 3.3 ks signal. This 
demonstrates that the XWT is picking out a well-defined, persistent signal, and will not 
misleadingly suggest a signal where there is none. 

We did not find any significant detections for the other nine AGN in our sample. No 
features had significances that exceeded the 99% confidence levels for both white noise and 
red noise processes (see Figs. [HI [ID]) and appeared at either too short (i.e. at timescales 
shorter than 3.6% the length of the observation) or too long (i.e. at timescales greater than 
half the length of the observation) a timescale. Some of the AGN in our sample have been 
studied before and previous reports of QPOs exist in the literature. We will discuss those 
results in more detail in Section 14.31 

3.3.2. Structure Function Results 

After correcting for pile-up, we subtract a constant level corresponding to Poisson pho- 
ton noise from the structure functions (Figures [T2 l[T3|) . The slopes are measured by fitting a 
power-law to the SF curve using the least-squares method in log{SF) —log{6t) space. Slopes 
are listed in Table [2] along with the characteristic time-scales of variability, which were mea- 
sured by identifying the times of turnover from plateau to power-law and vice versa in the 
SF curve. All of our structure functions have a flat plateau at short timescales corresponding 
to Poisson photon noise, most have a power-law portion, and some have a plateau at long 
timescales. We include light curves in Figures [HI [15] for comparison with the structure 
functions. 

The structure functions for all four observations of 3C 273 are shown in the first four 
panels of Figure [T2l The observation with the 3.3 ks quasi-period (upper left. Figure [T2|) is 
dominated by whitish noise around 3000s, as inferred from its flat slope; however, the SF 
is unsuited to quantifying the autocorrelation coefficient precisely. Recall that the wavelet 
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analysis finds an autocorrelation coefficient of a = 0.14, relatively small, and consistent with 
a fiattish structure function. We note that this observation also has the greatest excess of 
such noise above the photon noise, compared to the other three observations, consistent with 
this being a unique time series out of all those analyzed. 



Discussion 



4.1. Mass Estimates of 3C 273 



There are several mass estimates for 3C 273 obtained from different methods. One 
method is reverberation mapping whereby one uses the time lag of the emission-line light 
curve with respect to the continuum light curve to determine the light crossing size of the 
broad line region (BLR) and then assumes Keplerian condit ions in the broad line region gas 
motion (i.e. Mbh = v^Rblr/G) JPeterson fc Wandell[200o[ ). 



Reverberation mapping results based on the optical continuum (i.e. Balmer lines) place 



the mass of the central black hole in 3C 273 at 2.35toil ^& flKaspi et al.ll2000h . In a 



different study, iPian. Falamo. fc Treves! (120051 ) use Hubble Space Telescope UV luminosities 
to find the broad line region size. To do so, they derive a relationsh i p bet ween Rblr and 
UV luminosity using the er npirical relationship found bvlKaspi et al.l (120001 ) between Rblr 
and the optical luminosity. Plan. Falamo. fc Treves ( 20051 ) obtain a mass of 4.0^2 X IO^Mq 
for 3C 273, consi s tent with the Kaspi et al. (1200*0 ) value within errors. In another study, 
Paltani &: Tiirlerl (120051 ) look at the strongest broad emission UV lines (Lya and C IV ) in 



archival International Ultraviolet Explorer observations and obtain a mass of 6.59lg g^ x 
IO^Mq for the central supermassive black hole in 3C 273. 

There are also mass estimates for 3C 273 that do not come from reverberation mapping. 
Liang fc Liul (120031 ) find a black hole mass of 2 x 10^ Mq by generalizing the Elliot-Shapiro 
relation to the Klein-Nishina r egime for 3C 2 73 's g amma -ray flux obtained from EGRET. 
Another method is to use the iMcLure &: DunlopI (l200ll ) correlation between host galaxy 
luminosity and black hole mass wh ich obtains a mass of 1.6 x 10^ Mq with an uncertainty 
of 0.6 dex JWang. Luo. fc H0II2OO4I ). 



4.2. Underlying Physical Process for the QPO in 3C 273 



If the 3.3 ks quasi-period in 3C 273 represents an orbital timescale originating near a 
last stable orbit, it implies a central black hole mass of 7.3 x 10^ Mq for a non-rotating black 
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hole or 8.1 x lO^M!© for a maximally rotating black hole. These numbers agree with the 
Liang &: Lid (120031 ) mass estimate of 2 x 10^ Mq. However, these masses are substantially 
lower than those expected for supermassive black holes. 



The iPian. Falamo. fc Treves! (120051 ) estimate for the mass of the black hole in 3C 273 
at 4.0 X 10^ Mq points to an orbital period of ~ 200 k s for a last stable orbit of 3Rs and 
a period of ~16 ks for O.GRs for a rotating black hole. iPaltani fc Tiirlen (120051 ) estimate a 
mass for 3C 273 of 6.59 x 10^ Mq, which points to an orbital period of 3000 ks for a last 
stable orbit of 3Rs and a p eriod of 270 ks for O.QR s . The 3.3 ks quasi-period we find here 
is o nly about 2-20% of the iPian. Falamo. fc Treves! (12005! ) orbital timescale and 0.1-1% of 
the iPaltani &: Tiirler! (120051 ) orbital timescale, suggesting that this X-ray quasi-period is not 
caused by dynamical motion in the inner accretion disk. Furthermore, the inverse scaling 



between frequency and black hole mass yields an expected period o f t ~ SOOMbh/ 



based on the representative of HF QPOs in X RBs, GRO J1655-4 (lOrosz fc Bailvn 



Remi 



lard et al. 



1999 



Abramowicz et al.l |2004J ) . Using either the jPian. Falamo. fc Treves 



1997 



(12005! ) or IPaltani fc Tiirleii (120051 ) mass estimates yields a period that is one to two orders 



of magnitude higher than what we observe. 

Previous work has suggested that oscillations can occur in the innermost region of 
relativistic accretion disks due to their instability aga inst axisymmetric radial oscillations, 
possibly due to a magnetic field (IKato &: Fukud Il980l ). This has been proposed for X-ray 
binary systems, but the physical mechanis m responsib l e for these oscillations can be applied 
to other accretion disk systems like AGN. iPerez et al.! (Il997l ) analyze modes of oscillation in 
terms of perturbations of the general relativistic equations of motion of perfect fluids within 
the Kerr metric. They look at the case of a thin accretion disk around a Kerr black hole in 
order to determine black hole mass and angular momentum for different trapped modes. 

We propose that a g-mode oscillation of m > 3 is responsible for the 3.3 ks quasi-period 
in 3C 273. A g-mode (inertial) oscillation can be characterized as a restoring force that is 
dominated by the net gravitational-centrifugal force. These modes are the most relevant 
observationally since they appear to occupy the largest area of the disk and hence should be 
the most observable trapped modes (jPerez et al.lll997l ). 



Equation 5.4 of iPerez et al.l (119971 ) shows that the frequency of a quasi-period should 
be observed at 

/ = 7U{MQ/M)F{a)Hz (10) 

where a is the angular momentum parameter and M = M^gn- For m = and a = 0, 
F(0) = 1, while F{amax) = 3.443 where a max = 0.998 . This gives a mass that is too low. 
For m = 3, F{amax) ~ 59 (see F i gure 5 of iPerez et al.! (11997! )) and this gives a mass for 3C 



273 of 1.4 X IO^Mq. IPerez et al.l (11997! ) do not look at modes higher than 3 
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4.3. Previously Reported QPOs for AGN in Our Sample 



Fiore et al.l (119891 ) report a QPO in NGC 4151 around 5.8 k s with > 99% significance 
based on three EXOSAT observations. IVaughan &: Uttleyl (l2005bl ) reanalyzed these data sets 
and found that after fitting the red noise significance and Poisson photon noise components 
of the spectrum, the QPOs fall below the 95% threshold. Our XMM-Newton observations 
show a ~4.8 ks feature. This appears in our 57 ks observation (Obs. ID: 112830201) with 
96% red noise significance and in our 30 ks observation (Obs. ID: 112310101) with 99.4% 
red noise (see Fig. [TOl) . Even though this signal rises above the 99% red noise level in this 
observation, we discount it because it appears as part of a larger power structure in the 
GWPS and is not a well-defined peak. 



For NGC 5548, iPapadakis fc Lawrence! (119931) claim a 500s QPO in five out of eight 
EXOSAT observations. iTagliaferri et al.l (119961 ) reanalyzed the same data and found that one 
observation had detector problems. Also, in every case they found less than 95% significance 
by taking into account the uncertainties in modeling the spectrum. In our XMM-Newton 
data of NGC 5548, we report a 500s feature, but it has only a 93% red noise significance, 
and is seen in only one of our two observations (Obs. ID: 089960401, Fig. [TOj) . 



For MKN 766, iBoller et al.l (120011 ) claim a ~4200s QPO in a 30 ks XMM-Newton obser- 
vation. In our 128 ks observation, taken a year later, we see a signal at 4200s with 99.5% red 
noise significance, but it is dwarfed in the global wavelet power by a much stronger, wider 
broad peak (see Fig. [10]), possibly due to a secular change in flux over the observation. 



We do not de tect any signiflcan t feature for MCG-6-30-15 (ILee et al.ll2000f ). MKN 421 
or PKS 2155-304 (lOsone et al.l 1200 ll ). which have previously reported QPOs. To the best of 
our knowledge there are no published QPO claims for any of the other objects in our sample: 
3C 273, IRAS 13349+2438, NGC 3516. 



Halpern et al.l (120031 ) reported the discovery of a 2.08 day quasi-period in the NLSl 



galaxy T onS180 wi t h a 3 3 day observation taken with the Extreme Ultraviolet Explorer 
(EUVE) . IVaughanI (|2005l ) suggest that this periodogram is oversampled and so the signifl- 
cance is overestimated. Our wavelet analysis of this data shows a ~ 2 d period in the global 
wavelet power spectrum which rises above the 99% white noise level, but it has only 89.5% 
red noise signiflcance (Fig. [TUl) . Our wavelet analysis flnds an a of 0.8 for this observation 
and the structure function shows that red noise dominates, implying that this 2 d feature 
should be compared to red noise signiflcance and so is not signiflcant. 
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5. Summary &: Conclusions 

We applied the wavelet analysis technique to XMM-Newton observations of 10 AGN 
and detected a candidate 3.3 ks period in 3C 273. The cross wavelet transform shows that 
the 3.3 ks signal is present throughout the entire observation. 

If the 3.3 ks quasi-period in 3C 273 represents an orbital timescale originating near a 
last stable orbit, it implies a central black hole mass of at least 7.3 x 1 0^ Mq which does 
not agree with reverberation mapping mass estimates. Kaspi et al. ( 200o[) es timate the mass 



of the black hole in 3C 273 at 2.35 x 10^ Mq and IPaltani fc Tiirlei] (j2005|) find a mass of 
6.59 X 10^ M0. This suggests that this X -ray quasi-period is not caused by dynamical motion 
in the inner accretion disk. 

We suggest that oscillations with modes of three or highe r are occurring in t he accretion 



disk of 3C 273, producing the detected 3.3 ks quasi-period. iPerez et al.l (119971 ) shows that 
for m > 3 and maximum angular momentum one can obtain a mass for 3C 273 of 1.4 x 10^ 
Mq, consistent with the lower mass estimate obtained from reverberation mapping. 



We thank the anonymous referee for constructive comments that helped improve the 
paper. We thank Margo AUer, Robert Fender, Luis Ho, Jereon Homan, Jon Miller, and 
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Fig. 1.— Pile-up measurement from the SAS task EPATPLOT for MKN 421 Obs. ID: 
0099280101. The expected fractions of pixel events are the solid lines and the measured 
pixel events are the histograms (s,d,t,q stand for single, double, triple, and quadruple events 
respectively). In the range 0.75 to 10 keV we see that a larger than expected fraction 
of double events (third histogram from top, dark blue in electronic edition) is measured, 
indicating that pile-up is present in this source. [See the electronic edition of the Journal for 
a color version of this figure.] 
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Fig. 2. — Continuous wavelet transform for a sinusoid of varying frequency (top panel). 
Periodic behavior shows up as a pattern over all translations (y-axis) at a specific dilation 
(x-axis) in both the real part of the transform (middle panel) and the power (bottom panel). 
Here there are signals at 3s and 6s. Also given in the figure header are a (the autocorrelation 
of the signal), the range of dilations explored, and Ni (the number of dilations used). The 
hatched area represents the cone of influence. [See the electronic edition of the Journal for 
a color version of this figure.] 
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Fig. 3. — Global wavelet power spectrum for a sinusoid of varying frequency. The solid line 
is the power spectrum of the signal, which is compared to the power spectra of white and 
red noise random processes (broken lines). The 3 s and 6 s detections exceed the expected 
levels of white and red noise at the 99% significance level. 
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Fig. 4. — Cross-wavelet transform for a sinusoid. We use the cross-wavelet to compare a 
mock sinusoidal signal with a period of 6 s against the sinusoid with varying frequency from 
Fig. [21 The concentration in the real and power plots show that a 6s period exists in the first 
half of the data. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 5.— Effect of pile-up on structure functions of MKN 421 Obs. ID: 0099280101. The flat 
portion of the structure function curve affected by pile-up (left panel) falls below log{SF) — 
1, which corresponds to Poisson photon noise. In the right panel the central core of the 
source has been removed (i.e. pile-up is greatly reduced) and the SF curve is then at the 
expected value for Poisson photon noise (right panel). 
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Fig. 6. — Continuous wavelet transform for 3C 273 Obs. ID: 126700301. The real part of 
the transform shows an oscillatory behavior at 3.3 ks. There is also a concentration in the 
power plot at 3.3 ks. The 3.3 ks signal is circled in the both panels. The light curve (top 
panel) is binned up from 5s to 100s for clarity. The a for the unbinned 5s data is 0.14. The 
hatched area is the cone of influence: the region where edge effects become important [See 
the electronic edition of the Journal for a color version of this figure.] 




Fig. 7.— Global wavelet power spectrum for 3C 273 Obs. ID: 126700301. The solid line is 
the power spectrum of the signal, which is compared to the power spectra of white and red 
noise random processes (broken lines). The 3.3 ks detection exceeds the expected level of 
red noise with a 99.979% probability of detection. 
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Fig. 8. — Cross-wavelet transform for 3C 273 Obs. ID: 126700301. We use the cross-wavelet 
to compare a mock sinusoidal signal with a period of 3282 s with the observation. The 
concentration in the real and power plots show that a ~3300s period exists in the data. [See 
the electronic edition of the Journal for a color version of this figure.] 
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Fig. 9. — Global wavelet power spectra for all observations in our sample with log(Global 
Power Spectrum). The period is in seconds. 



-27- 




Period 



Fig. 10. — Global wavelet power spectra for all observations in our sample with log(Global 
Power Spectrum). The period is in seconds except for TonSlSO where the period is in days. 
TonSlSO is not in our sample, but is discussed in Section 4 and so we include it here. 
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Fig. 11. — Histograms from Monte Carlo significance tests. The histogram shows that for 
simulated light curves created with Poisson photon noise there are four false detections (i.e. 
periodicities at a significance of > 99.99%) at timescales 2000 < r < 1.5 x lO^s, translating 
to a probabihty of 0.4%. The y-axis is cutoff at 45 but rises to 650. 



-29- 



I I 

cn 

O 



rrrq i i i i i iiij i i i i i iiij i i i i i iiij 

'^■ife- 3C273 Oba.lDrOI 26700301 



0,1 



u 
c 



I II I ! I I I I I 




mr] I I I I I iiij I I I M I iij I I I I I iiij I r 

- 30273 Obs.lDll 26700001 




Tmj 1 — I I I I iii| 1 — I I 1 1 iii| 1 — I I 1 1 in| 

r- 3C273 ObsJD: 136550101 



LLUl I I I I I - I I 



TO fSe uJnb Tor 
rTTT| 1 — I I I I iiij 1 — I I 1 1 1 iij 1 — I I I 1 1 iij r 

IRAS 13349+2433 
Ob9,ID096010101 




0,01 



TTTTj ! — I M I Ml| 1 — I M I Mlj 1 — I M I I! I| 1 — n 

3C273 Obs,ID;01599S0101 



LLul I 1 I -J I |. | lliL ■■ .■:1?>>J-'iiilrJ- 



10 



WD 



1000 




10* 




:0,1 



T| 1 — I 1 1 II iij 1 — I I I iiiij 1 — I I I iiiij 

WS1 Obs.lD:1 1 IBDOIOf 



" ' 4 ' " ' 



"4" ' ' -'^ ' ^i^'^ 




0.1 =- 



□.[}i 



^ Mlj 1 1 1 1 1 Mlj 1 1 1 1 II M| 1 


1 




^ IVICG-6-30-15 






! 0bsJD:02974Q7(5l 












1 'III < ■ 1 1 1 


1 1 


1 1 1 It- 



10 



1Q0 



IDOQ 



IC 



■|gpnr| 1 I I I I Mlj 1 I I 1 1 Mlj 1 — I I 1 1 Mlj r 

MCG-6-30-15 
ObsJD:029740f 



1.1 




lOr 



TTTT| 1 I 1 1 I I III l| r 

MKN 421 



0,1 
D.Ql 



r Qbs.lD:0992SOT01 



iB^ 




-1 1 lll 1 1 1 1 1 1 Illj 1 


1 1 1 1 1 Ml 




r WKN 421 






= 0bs.lD:0992a0301 






' ' iil i I I' 1 " 1 1 


1 





TO* 

Log [Time(s)] 



Fig. 12. — Structure functions for each object in our sample. Poisson photon noise has 
been subtracted. All of our structure functions have a flat plateau at short timescales 
corresponding to Poisson photon noise and most have a power-law portion corresponding to 
red noise. 
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Fig. 13. — Structure functions for each object in our sample. Poisson photon noise has 
been subtracted. All of our structure functions have a flat plateau at short timescales 
corresponding to Poisson photon noise and most have a power-law portion corresponding to 
red noise. 
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Fig. 14. — Lightcurves for the 10 AGN in our sample. All observations are in the energy 
range 0.75 to 10 keV. Most observations are binned to 100s except for MKN 421, NGC 5548 
(Obs ID: 089960401) and IRAS 13349+2438 which are binned to 50s. 
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Fig. 15. — Light curves for the 10 AGN in our sample. All observations are in the energy 
range 0.75 to 10 keV. Most observations are binned to 100s except for MKN 421, NGC 5548 
(Obs ID: 089960401) and IRAS 13349+2438 which are binned to 50s. 
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Tablc 1. Log of Observations 



Object 


Observation ID 


Start Date 


Length (ks) 


(Oounts) 


30 373 


126700301 


2000-06-13 


66 


140 


30 373 


126700801 


2000-06-17 


60.6 


120 


30 373 


136550101 


2003-01-05 


88.6 


160 


30 273 


159960101 


2003-07-07 


58 


230 


IRAS 13349+2438 


096010101 


2000-06-20 


44.6 


10 


M81 


111800101 


2001-04-22 


130 


20 


MOG-6-30-15 


029740701 


2001-08-02 


127 


70 


MOG-6-30-15 


029740801 


2001-08-04 


125 


120 


MKN 421 


099280101 


2000-05-25 


32.5 


740 


MKN 421 


099280301 


2000-11-13 


46.6 


1060 


MKN 766 


109141301 


2001-05-20 


128.5 


90 


NGO 3516 


107460601 


2001-04-10 


129 


20 


NGO 3516 


107460701 


2001-11-09 


128 


12 


NGO 4151 


112310101 


2000-12-21 


30 


20 


NGO 4151 


112830201 


2000-12-22 


57 


25 


NGO 5548 


089960301 


2001-07-09 


93.4 


75 


NGO 5548 


089960401 


2001-07-12 


37 


90 


PKS 2155-304 


124930201 


2000-05-31 


59 


280 


PKS 2155-304 


124930301 


2001-11-30 


44.6 


380 



Note. — 



Events are grouped in 5s bins. 
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Table 2. Structure Function Slopes and Timescales 



Object 


Observation ID 


Slopes 


Transition Time (s) 


Turnover Time (s) 


30 373 


126700301 


1.23 


1.5 X 10* 




30 373 


126700801 


2.09 


1.5 X 10" 




30 373 


136550101 








30 273 


159960101 


1.68 


2 X 10* 




IRAS 13349+2438 


096010101 


1.68 


1000 




M81 


111800101 


0.95 


7000 




MCG-6-30-15 


029740701 


0.82 




10" 


MOG-6-30-15 


029740801 


1.11 


200 




MKN 421 


099280101 


1.25 


400 




MKN 421 


099280301 


1.13 


350 


7000 


MKN 766 


109141301 


0.67 


100 


3 X 10"' 


NGC 3516 


107460601 


1.18,1.39 


2000 




NGO 3516 


107460701 


1.86 


2 X 10'' 




NGO 4151 


112310101 








NGO 4151 


112830201 


1.19 


7000 




NGO 5548 


089960301 


0.99,1.64 


1200 




NGO 5548 


089960401 


0.097,2.75 


1000 




PKS 2155-304 


124930201 


.99,.59,.45,1.69 


2000 




PKS 2155-304 


124930301 


1.59 


1100 





Note. — Transition time corresponds to the time at which the SF curve changes from plateau to 
power-law. Turnover time is the time at which the power-law portion of the SF curve changes to a 
plateau. 



